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Abstract 

In this paper we find the most general self-similar, homogeneous and isotropic, Ricci flat cosmologies in 5D. 
These cosmologies show a number of interesting features: (i) the field equations allow a complete integration in 
terms of one arbitrary function of the similarity variable, and a free parameter; (ii) the three-dimensional spatial 
surfaces are flat; (iii) the extra dimension is spacelike; (iv) the general solution is Riemann-flat in 5D but curved 
in 4D, which means that an observer confined to 4D spacetime can relate this curvature to the presence of matter, 
as determined by the Einstein equations in 4D. We show that these cosmologies can be interpreted, or used, as 
5D Riemann-flat embeddings for spatially-flat FRW cosmologies in 4D. In this interpretation our universe arises 
as a topological separation from an empty 5D Minkowski space, as envisioned by Zeldovich. 
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1 Introduction 



The question of whether there was a Beginning of the Universe, is a truly challenging puzzle for physics and for 
philosophy 1 . In classical four-dimensional general relativity, the singularity theorems imply that the big bang was 
a unique birth event for a universe filled with matter satisfying the weak, dominant and strong energy conditions 
2 - 3 . The big bang model postulates that our observable universe originated in a singularity sometime between 10 
and 20 billion years ago. Before the bang nothing existed, not space, time, matter, or energy. 

Astronomical observations support the notion that our universe began at some point in the finite past 0]-[B]. 
The question is how did it begin. Twenty seven years ago, in a paper concerning the idea of a spontaneous birth 
of our universe, Ya. B. Zeldovich observed that "there is a certain arbitrariness and fuzziness in the very concept 
of spontaneous birth" [7] . He inquired whether spontaneous birth emerges (i) "out of nothing" , or (ii) in a space of 
more dimensions, or (iii) as a topological separation from an initially given empty Minkowski space. 

Today, in the literature we find a huge number of papers examining theories for the birth of our universe along 
the lines (i) and (ii) envisioned by Zeldovich3. In fact, there are so many interesting papers and books about these 
topics that for us here is impossible, and far beyond the scope of this work, to give a thorough list of references. 
Therefore we restrict ourselves to mention just few representative works where one can find a detailed bibliography: 
tunneling from nothing [11]- [17]; spontaneous creation of the braneworld [H]-[2T]; induced matter in Kaluza-Klein 
gravity and STM (Space-Time-Matter) theories [22] -[25]. 

Concerning the third alternative mentioned by Zeldovich, the hypothesis that singularities in AD are induced 
by the separation of spacetime from the other dimensions has been examined by the present author [27] ; Seahra 
and Wesson have thoroughly investigated the structure of the big bang from a five-dimensional embedding [55] for 
the standard spatially- flat AD FRW models [28]. Besides these studies, there are only few more works related to 
this alternative [29]- [33]. The aim of this paper is to present a general class of spatially homogeneous and isotropic 
solutions of the Einstein field equations in 5D that is relevant to this alternative. 

To be more precise, in many Kaluza-Klein and braneworld theories the higher dimensional space is assumed to 
be either Ricci-flat or anti-de Sitter. In this work we develop a family of metrics that are Ricci-flat and Riemann-flat 
in 5D, but whose four-dimensional subspaces are curved. Therefore, they are all equivalent to an empty Minkowski 
space in 5-D. However, since the Riemann tensor of four-dimensional subspaces is non-vanishing, for an observer 
confined to 4-D the spacetime is not empty but contains matter as determined by the Einstein equations in AD. 

The solutions arise from the observation that cosmological models are self-similar [34] . In order to illustrate this 
in 5-D, let us consider the metri (3 [22] 

dS 2 = Cy 2 dt 2 - Dt 2 ' a y 2 l^ \dr 2 + r 2 (d0 2 + sin 2 6d0 2 )} - , CaH * dy 2 , (1) 

(1 — a) 

where a is a free parameter and C as well as D are constants with the appropriate units. The properties of this 
metric have widely been discussed in the literature, e.g., [24) . [28] . [35] . It is Riemann-flat and apparently empty in 
5-D, but reduces to the well-known AD FRW models with flat 31? sections on y = constant hypersurfaces (henceforth 
denoted by S y ). It is easy to verify that {J) admits a homothetic Killing vector, viz., 

C (9AB = 2 9AB , with ? = (^L_ y r, 0, 0, g^>) . (2) 

With the transformation of coordinates 

t _ f a/(2a-l) j r ^ f> y _ >g (a-l)/(2a-l) j (3) 

1 It should be mentioned that in addition to the theories envisioned by Zeldovich, there are many other ones, e.g., evolution of our 
universe from a 4D Minkowski spacetime [8], [9]. More recently, we find black hole generation of universes; universes spontaneously 
creating other universes; creation of universes by intelligent life. For a wonderful popular review, see for example |10| 

2 Conventions: Throughout the paper we use geometric units where c = G = 1; t = x°, r = x 1 , 9 = x 2 and <j> = x s are the usual 
coordinates for a spacetime with spherically symmetric spatial sections; y = x 4 represents the coordinate along the extra dimension; the 
signature of the 5-D metric is (+, — , — , — , e) where e can be either —1 or +1 depending on whether the extra dimension is spacelike or 
timelike. The range of tensor indices is A, B... = — 4 and fj,, v, ... = — 3. 
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we have 

C A -> Z A = (i, f, 0, 0, y) (4) 

and the metric |T]) becomes 

dS 2 = A ^(l-a)/(2a-l)^"2 _ g f/(2a-l) [ df 2 + r 2 ^2 + ^2 ^2)] _ ^ ^ Q /(2 Q -1) ^-2 _ (g) 

where A and £> arc dimensionlcss constants and 



The above clearly illustrates the self-similar nature of (JTJ) . In this work we extend our previous studies [32] and 
[54] to obtain a general class of self-similar cosmologies that are Ricci-fiat and Riemann-flat in 5D, but curved in 
AD. Consequently, they can be used, or interpreted, as 5D Riemann-flat embeddings for our AD universe, which is 
consistent with Zeldovich's notion that our universe could have emerged as a topological separation from an initially 
given empty Minkowski space in more than four dimensions. 

The paper is organized as follows. In Section 2, we deduce the equations with self-similarity in 51?. In Section 
3, we obtain the general solutions to these equations; they contain an arbitrary function of the similarity variable, 
and a free parameter. In Section 4, we focus our attention on a class of solutions that admits a particularly simple 
homothetic Killing vector in 5D, and study their possible application in AD. We use the flexibility of the solution 
to develop some cosmological models in AD. In Section 5, we present a summary of our results and propose several 
extensions for this work. Finally, in the Appendix, we discuss the homothetic symmetry on E y . 



2 Field equations 

In this section we present the field equations for self-similar cosmologies in 5D. We start with the line element for a 
spacetime that has spatial spherical symmetry 

dS 2 = e^'^dt 2 - e Mr ^' y) dr 2 - R 2 (r, t, y) [d6 2 + sin 2 9d<j) 2 ] + e$ 2 (r, i, y)dy 2 . (7) 

For cosmological models we assume spatial homogeneity, which means that the metric is invariant under spatial 
translations. However, we do not make any assumption regarding the curvature of 3-D-space. Now, as illustrated 
above ([5]) , by a suitable transformation of coordinates in a self-similar model all the dimcnsionless quantities can be 
put in a form where they are functions only of a single variable (say £) [36]- [43]. Thus, in the case under consideration, 
in "self-similar" coordinates i, f, and y, the line element ([7]) can be written as 

dS 2 = e u ^dP - e^df 2 - fV l( ^ [d6 2 + sin 2 9d<j) 2 ] + e$ 2 (£)dy 2 . (8) 

On the other hand, we have no reason to set £ = t/y as in ©. In principle £ can be any function of t and y, namely 

£ = (9) 

With this choice the line element ([8]) is self-similar but not necessarily admits a homothetic Killing vector (More 
comments about this at the end of Section 3). In what follows we are going to suppress the bar over the self-similar 
coordinates. 

The metric functions in ([7]) and ([5]) have to satisfy the 5D Einstein field equations in apparent vacuum, which in 
terms of the Ricci tensor are 

Rab = 0. (10) 

At once we note that R i = required 

A € =/x 4 (11) 

from which we get 

e M«) =C e^\ (12) 



3 In what follows fp denotes derivative of / with respect to 5; dots and primes stand for derivatives with respect to t and y, respectively. 



3 



where Cq is a constant of integration. Now, from R\ = R\ = i?| it follows that 

Co = 1. (13) 

Consequently, the requested self-similarity 1(5)), together with the field equations (fit))) , demand i? = re A / 2 . Therefore, 
the line element for self-similar cosmological models in 5D reduces to 

dS 2 = e v ^dt 2 - e A(5) [dr 2 + r 2 (d9 2 + sin 2 9d4> 2 )} + e$ 2 (£)dy 2 . (14) 

This implies that in self-similar cosmologies the 3D spatial sections (t — constant, y — constant) are flat, in agreement 
with astrophysical data [1]-[B]. For this metric the Ricci tensor in 5D has four non-trivial, independent, components. 
These are 

Ron, R04, Rn = — 5- = „ . n ~ ) -^44- (15) 
r z H sin 

A simple analysis of equation R 04 = 0, indicates that self-similarity requires that the ratio [£'/(£'£)] be some function 
of £. Clearly, any separable function of t and y will do the job. Therefore, without loss of generality we can set 

t T(t) 

^WY (16) 

where T and Y are some functions to be determined by the field equations in 5D. Consequently, we have four 
differential equations for five unknown, viz., 

A(0, T(t), Y(y). (17) 

We note that e is also to be defined from the field equations. We will see that we can solve the field equations (|10p 
in terms of one arbitrary function and a free parameter. The four equations to be integrated are: 



1. -Roo = 



e $ 2 T 2 (<*« + 3A 2 - a*A« + - + 2e$ 2 f y ( 3 A, + *§* 

y'Vf (2^ + ,f + 3^A S - ^i) + 2 y'V (2 - ^) to = 0, (is) 



2. i? n4 = 0, 



2A« + A 2 + ^i-^-^i=0. (19) 



3. Jin = i? 22 /r 2 = R^/{r 2 sin 2 0) = 0, 



e$ 2 T 2 ^2A« + 3A| - ^A 4 + + 2e$ 2 f yA e + 

y' 2 e^ 2 (2A« + 3A| + ^ - ^) + 2y'V (2 - ) ^ = 0. (20) 



4. Finally, i? 44 = 0, 



2e$ 2 T 2 f 2*« _ + 3>£A + 4e$2f y (*£ , + 

v $ $ $ / V * 



V V I2 C 2 fo,, 1 ,,2 2f£$£ 2 



e'Y' 2 e [ 2v a + v{ + 6A« + 3A 2 - I + 2y' 2 e " I 2 - ) £ fa + 3A C ) = 0. (21) 
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3 Integrating the field equations 

Let us first notice that (fT§|) can be easily integrated as 

£ 2 A 2 = c 2 $V"- A \ 



(22) 



where C is a constant of integration. On the other hand, the field equations (fT5)) . (|2T)|) and (f2"Tj) have the following 
structure 



£s) m + f P 2 " ) G(0 + ^ + 7(0 ( 2 " W 1 = 11 



(23) 



where P 1 , G, i? and 7 symbolize the corresponding functions of £ in these equations. Therefore, in order to preserve 
the self-similar symmetry, we have to require 

YY" 

2-^ = 1, (24) 



y/2 

where I is a separation constant. Integrating this expression we obtain 

Y' = aY^ 2 - l \ 
where a is a constant of integration. Thus, we find 

/ ty\ 



yi2 



)- 


f 2 




a 2 T (4-2l) 



^- 2l \ and 



Y'2 



y 2 T (3-2l) 



(3-2/) 



(25) 



(26) 



Consistency of (I23[) demands the quantities inside the square brackets to be constants, which requires 

tV('-i) for I £ 1, 

e 13 * for I = 1, 

where (3 is some constant. Consequently, without loss of generality we can set 



T 



(27) 



i = 



*) for I ^ 1, 



(28) 



for I = 1. 



Let us now calculate (Rqo — -R44), and use (fT9|) to express the second derivative A{{ in terms of the first derivatives 
of v and $. As a result we obtain 

(29) 



( TT — T 2 ) = (Y rz - YY") 



12 \r\rll\ c2v 



Solution for 1^1: In this case (|2T)|) yields 



$ 2 = (~e)e 



(30) 



From this expression, it follows that the extra dimension should be spacelike, i.e., e = — 1. Now, from (|22p and (f3T)|) 
we get 

= (l\^-i) Xi e x ' 2 . (31) 



In summary, we have found that the line element (|14| with 



, and e = — 1, 



(32) 



is a solution of the bD Ricci-flat equations for any value of the parameter I ^ I, constant C, and arbitrary function 
A = A(£). In addition, it admits a homothetic Killing vector in 5-D, viz., 

C^g AB = 2g AB , with £ A = (t,r,0,0,y). (33) 

It should be noted that (j3"2")l includes the family of 5D metrics given by |T]), or © for the particular choice e x ^ = 

B £2/(2 a -i) and l = 2 

Solution for 1 = 1: In this case the field equations require j3 = a. Therefore, the solution is 

i)eV^/2 ; <j> 2 = e K« e=(^), and e = -l. (34) 
It is obvious that this solution does not admit a simple 5D homothetic vector as ([2]) or (|33|) (See bellow). 
3.1 The Riemann tensor in 5D 

It can be verified that, for both solutions, all the components of the 51? Riemann tensor Rabcd vanish identically. 
Therefore, they are equivalent to an empty Minkowski space in 5D (A4§). Consequently, there exist some coordinate 
transformation 

T = r(t,r,y), R = R(t,r,y), ip = i>(t,r,y), (35) 
that brings (f3"2")l and (f3"4")) to the line element in 

dS 2 = i] AB dx A dx B = dr 2 - dR 2 - R 2 (dO 2 + sin 2 9d<t> 2 ) - d^ 2 , (36) 

in Minkowski coordinates x A — (r, R, 9, 0, ip). For the particular solution ([4]), the explicit coordinate transforma- 
tion from ([I]) to (|36|) is known and has been amply discussed in the literature, see for example [28] and references 
therein. 

Before going on, and in order to avoid misunderstandings, we should comment about the homothetic nature of 
the above solutions. Certainly, any 5D vector of the form 

tfi = {[r + h{RM> R, 0, 0, [tl> + f 2 {T,R)]}, (37) 

where f\ and f2 are arbitrary functions of their arguments, is a homothetic Killing vector for (|36|). viz., C^ M rjAB = 
2i]ab- Thus, the Riemann-flat solutions (132|) and (|34|) admit an infinite number of homothetic killing vectors, which 
correspond to the nondenumerable infinity of choices of /i and /2 in (|37|) . 

In this context, the difference between (|3"2"|) and ([3^)1 is that, among all the possible choices, in self-similar coordi- 
nates t, r and y there exists a very simple homothetic killing vector for the first solution , namely f = (i, r, 0, 0, y), 
while for the second solution the homothetic vectors (in self-similar coordinates) are much more complicated. 

4 Interpretation in AD 

In this section we discuss possible four-dimensional interpretations of the 51? line element (132[) . For simplicity we set 

I = 2, (38) 

in such a way that the similarity variable takes the form £ = t/y, as in ([6]). Thus, in what follows we consider the 
5-D metric 

dS 2 = (L\ A ? e A <^ 2 dt 2 - e A «) [dr 2 + r 2 (d9 2 + sin 2 9dcf> 2 )] - (J^J £ 2 A c e A «>/ 2 dy 2 , (39) 
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with 

t-1. (40) 

The first important question is how to recover our AD spacetime from 5D. The most popular approach is to 
assume that our AD spacetime is a hypersurface S y : y = yo = constant, which is orthogonal to the 5D unit-vector 

along the extra dimension, although a dynamical foliation is also possible |44j . [45] . The second important question 
is how to construct the metric of the physical spacetime from the one induced on T, y . There are various approaches 
in the literature: (i) the canonical metric, which assumes $ = 1 and factorizes the AD part of the 5D metric by an 
y 2 term [46]- [49]; (ii) the conformal approach, where the AD part of the 5D metric is factorized by an $> N term (see 
for example |50] and references therein); and (iii) the one where the spacetime metric is identified with the metric 
induced on Sy. 

An exhaustive treatment of all possibilities is beyond the scope of this work. We present here an introductory 
analysis where we follow the approach (iii) mentioned above. The induced metric h a p on hypersurfaces T, y is just 
the AD part of the 5D metric (|59"1) . viz., 

ds 2 = h^dxPdx" = (J^j X^^dt 2 - e x ^ [dr 2 + r 2 (d9 2 + sin 2 <f>)] . (42) 

We immediately notice two things: 

Firstly, that the spacetime metric h a p is not self-similar along = (t, r, 0,0), which is the projection of £ A = 
(t, r, 0, 0, y) on E y , i.e., 

C u V ? 2V. fOT £p = (*> r, 0, 0) . (43) 
In fact, we show in the Appendix that there is only one family of homothetic solutions on E y . But the homothctic 
vector in AD, say is not parallel to the projected 

Secondly, although (|39[) is Riemann-flat, the hypersurfaces H y are curved. Indeed, the non- vanishing components 
of the Riemann tensor calculated on T, y are 



# 



oioi — 



#0202 


#0303 




r 2 sin 2 9 


#1313 


#2323 


sin 2 e 


r 2 sin 2 9 



(2A« + A 2 ) 



8y 



2 



#1212 = ^f = -r?T7=- \ 7 ■ ( 44 ) 
sir 9 r 2 sir 9 Ay 2 

Now, the requirement /loo 7^ implies 7^ 0. Therefore, R a pfw ^ 0. What this means is that an observer, who 
is confined to making physical measurements in our ordinary spacetime, can explain the curvature of T, y as being 

produced by ("effective") matter whose energy- momentum-tensor T^J is given by the Einstein equations in AD, 



where = (dg^ v / dy) /2. For the case under consideration, a simple calculation yieLdfl 



(45) 



tepWf) = ^ A 5 e- A «)/ 2 , (46) 
8np^ = U + ^) e-^/ 2 . (47) 



We now proceed to illustrate the above discussion with some examples. 



4 Here 'f> = h°°T^ ff) , and p^fD = -h 11 T^ 1 ff) = -h 22 T^ ff) = -h 33 T^ ff) . 
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4.1 Solution with $ = 1: a Riemann-flat embedding for the de Sitter universe 

Many authors use the five available degrees of coordinate freedom to set g^ = and $ = 1. This is the so-called 
"Gaussian normal coordinates system" based on Y, y , where n A is taken to be geodesic in 5D. In this system we can 
easily integrate (|32j) to obtain (with k = 2) 

1 .o C x2 



dS 2 = ±dt 2 -[B+j) [dr 2 + rW] - dy 2 (48) 



where B and C are dimensionless constants. With the transformation of coordinate^ 



-=-^, t^e- l l\ (49) 



t L 

where L is some constant length, (|48p becomes 

dS 2 = ^dP -(B + Cye 1 /^ 2 [dr 2 + r W] - dy 2 : (50) 

where the new constant C has dimensions of (length) -1 . On E y the effective matter is given by 

B (rfri 3C 2 e 2i / L B fpfn Ce i / L (2B + 3Cye i / L ) 

8np {eff > = = — , 87rn (e//) = i = — ^ ^. 51 

(B + Cye*/ L ) 3 yiB + Cyet/L) 2 

We note that 

87rp (e//) _> _8vrp (e//) = 4 as f oo. (52) 

Therefore PS]) , or (|50p , for t w is a Riemann-flat embedding for the de Sitter universe with cosmological "constant" 
A = 'i/y 2 . 

4.2 Solution with e v = 1: a Riemann-flat embedding for Milne's universe 

The choice /loo = 1 (and hoj — 0) is usual in cosmology; it corresponds to the so-called synchronous reference system 
where the time coordinate t is the proper time at each point. In such a system the line element (|39|) becomes 



dS 2 = dt 2 -{B + C£) 2 [dr 2 + r 2 dQ 2 ] - £ 2 dy 2 . (53) 
On Y. y the effective matter quantities are given by 

We note that the equation of state p = — 3p appears in different contexts: in Milne's universe; in discussions of 
premature recollapse problem [55]; in coasting cosmologies [S3]; in cosmic strings [S3], [5SJ; limiting configurations 
[5B] ; the exterior spacetime for stellar models in 5D Kaluza-Klein gravity [57] . 



5 Setting B = 0, and making a coordinate transformation dt/t = — yj A(t) /3 dt, where the function A(f) has units of (length) 2 , @8) 
can be written as dS 2 = y 2 Mp-dP — C 2 y 2 e 2 1 V A (')/ 3d ' j^ r 2 _|_ r 2^Q2j _ ^2^ wn ; cn j s identical to the 5D line element discussed by Bellini 
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4.3 Riemann-flat embedding for a AD universe filled with ordinary matter and a 
cosmological constant 

In order to obtain an equation for the unknown function A(£), let us assume that the effective matter can be separated 
as 

8irp eff = 87175 + A, 

87rp e// = Hnp-A, (55) 

where p and p are the density and pressure of the cosmological fluid and A is the cosmological constant. In addition, 
as it is usual in cosmology, we assume that the density and pressure satisfy the barotropic equation of state 

P = n P, (56) 

where n is some constant commonly restricted by \n\ < 1, which follows from the dominant energy condition [2], [3]. 
Thus, 

p eff +p eff A(n + 1) . eff eff . 

(n + 1) 8tt 



Using ([46]) and gT]) we find 



2A ee + (3n + 2)A| - /i A ? e A/2 = 0, with h = 4A ( n + 



(58) 



A= 77Z ?V 2 - (59) 



where h is a constant. Thus, in general 

h C 
4(n + l)y 2 

Setting 5(0 = e x ^/ 2 , ([53]) beco mcs 

2S% + 2(3n + l)S , f-/i S 2 ,% =0, (60) 

whose first integral is 

b 2Ay 2 

S S = + g(3 «+l) ' Where a =~3C~' (61) 

i.e., /lo = 6a(n + 1) and 6 is a constant of integration. Thus, in the case under consideration the spacetime part of 
the 5D metric (f3T)|) is given by 

ds 2 = 1 US 2 + s(1 ^ 3w) ) * 2 ~ 5 2 [dr 2 + r 2 (dfl 2 + sin 2 0# 2 )] . (62) 
The corresponding matter quantities arc 

3C6 , 3Ca 

P = 2y 2 ^(n+D ' P = UP > A= 2^- (63) 



Let us consider the deceleration parameter q, which is defined as 



(64) 



where S T = (l/y)S^e u l 2 represents the derivative of the scale factor with respect to the universal time r, which is 
related to the coordinate time t by the expression dr = e v l 2 dt. A simple calculation gives, 

_ in + 1 3a (n + 1) 

q ~~ 2[a + b/S^ n + 1 )]' { ' 
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In order to simplify this expression let us calculate the density parameters fl m = 8np/3H 2 and J7a = A/3H 2 , where 
H = S T /S is the Hubble parameter. We obtain 



Ca 



(b/a) 
53(1+1) 



fir. 



(6/0) 



Consequently, 



53(n+l) +(6/ )' 



SI, 



g3(n+l 



53(n+l) + (5/ a )- 



Substituting this expression into (|65p . and using that f2 m + f^A = 1, we get 



(66) 
(67) 



Q = 2^™(" + 1) - 1- 



We note that q changes sign at fi m = n£ rit) = 2fi A /(3n + l): for fi m < Q [ ,n"' the expansion is slowing down (q > 0), 
for f2 m > flm"^ the expansion is speeding up (q < 0). Setting n = 0, in concordance with the fact that our present 



(crit) 



(68) 



universe is matter-dominated (p = 0), and fi„ 
parameter today, viz., 

which is consistent with observations [3], [BJ. 



Q 



0.3 we obtain the approximate value of deceleration 



Q\today ~ —0.55, 



(69) 



4.3.1 Analysis of the evolution equation (|61 1> 

Let us notice that, in general, (|6ip cannot be integrated in terms of elementary functions. However the approximate 
behavior of the solution is as follows. In the very early universe, when 5 ~ 0, the second term in the r.h.s. of (|6ip 
dominates over the first one. Thus, we find S ~ ^ 1 /( 3 ™+ 2 ) which implies that the actual behavior of the scale factor 
depends on the equation of state. 



False vacuum: n = - 1. Following Zeldovich [7], one may imagine that after the spontaneous birth (where the 
spacetime separates from the extra dimension) our universe enters a de Sitter phase of expansion. In this phase 
n = — 1 and S ~ 1/t on every hypersurface E y , which means that S = at coordinate time t = 00. With the 
transformation of coordinates t — > e~ UT ' , where a; is a constant, we obtain S ~ e" 1 ". Thus, S = at r = — 00 which 
corresponds to t — 00. Thus, during the inflationary period the expansion of the universe is described by 

ds 2 = dr 2 - Ae 2 V A <^n/3r ^2 + ^ (70) 

where A( in fi) = 3oj 2 is the "effective" cosmological constant in this epoch. According to (|63|) it is related to the 
present cosmological constant by 

SirpWV = A (mfl) =A(l + ^. (71) 

However, it should be noted that the de Sitter solution is unstable under small perturbations. Therefore, it is 
impossible to extrapolate it to r = —00 7 . What this means is that (|70|) does not describe 5 = 0, which in singular 
on £ y (but non-singular in 5D), corresponding to the moment of spontaneous birth of our universal 



FRW evolution. The switch from the de Sitter exponential expansion to the radiation dominated FRW universe 
is a jump of pressure from p = —p to p = p/3. In this case, 5 ~ i 1 / 3 on every Ej,. Now changing t — > r 3 / 2 the 
scale factor becomes 5 ~ r 1 / 2 , which is the usual expression for a radiation dominated era. For any n / -1 the 
transformation t — > r 2 ( 3 ™+ 2 )/[ 3 ( 1 +")] allows us to recover the familiar spatially-flat FRW models, viz., 

ds 2 = dr 2 - Br 4 /[ 3 ( 1+ ")1 [dr 2 + r 2 dtf] . (72) 
6 Models for inflation estimate that the de Sitter phase of exponential expansions begins at r = 10 — 42 s 1581 . 
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Let us call To the moment of transition from the de Sitter phase to the relativistic plasma (n = 1/3) FRW solution. 
A jump in the pressure indicates that the second derivative of the scale factor is discontinuous, which in turn implies 
a jump in the deceleration parameter q; in the present case from q = — 1 to q = 1. However, in order to guarantee 
the continuity of energy density and the Hubble term, the scale factor and its first derivative must be continuous 
during the jumrj^- This allows us to relate tq and A( in y;) , namely, 



T0 = Wa • (73) 

1 V A o/o 



Our toy model does not provide any mechanism for calculating the time of transition, neither is the object of the 
present work, however various inflationary models suggest that inflation ends at To = 10~ 32±6 s [58]. Using this 
number as a reference we obtain 

A (in/0 ~ 10 48±12 cm- 2 , (74) 

which is huge compared to the present value of the cosmological constant A ~ 10~ 51 cm~ 2 , measured in the current 
ACDM model of cosmology Thus, 

A (m/() « 10 99±12 A. (75) 

Accelerated expansion. As the universe expands the deceleration parameter q, which is given by (|65|) . changes 
sign again at 



■5(5=0) - 



3n + 1 / \infi) _ 1 
A 



l/3(n+l) 

(76) 



where we have used (|T1 [) to express the ratio (b/a). This expression uncovers two interesting physical concepts: (i) 
if h-Unf) — A, then the universe would have never changed from decelerated to accelerated expansion, and (ii) the 
very large ratio A(j n m//A, evaluated in (I75j) . explains the very large radius of the universe. 

The present value of q, as calculated in (|69p. is approximately —0.55 which means that our universe is in a phase 
of accelerated expansion. At late times, for large values of S, (S > S/ q —Q\), the first term in (pTlj) starts dominating 
over the second one. Asymptotically, for S S> Sf g —Q\, in terms of universal time t our world enters a new era of 
exponential expansion with q = — 1 and 

ds 2 = dT 2 _ e 2^r ^2 + r 2 d ^ (?7) 

In summary, the toy cosmological models considered in this section allow us to conclude that that the Riemann-flat 
embeddings given by (|32[) are rich enough as to accommodate the essential features of the evolution of our universe. 



5 Concluding remarks 

In this paper we have found the most general self-similar, Ricci-flat, homogeneous and isotropic cosmologies in 5D. 
Self-similarity requires that all the dimcnsionless quantities in the theory be functions of a single variable £, which in 
the case studied here is some combination of t and y ©. As a consequence of the symmetry, and the field equations 
Rab = 0, we have found that (i) the three-dimensional spatial surfaces defined by t — constant and y — constant 
have to be flat (fT4|) ; (ii) the self-similar variable can be taken either as £ = t/y ([32]) . which is the more general case, 
or as £ = e at /e ay (|34|) : (iii) the extra dimension must be spacelike (|30|) . Then, we obtained the general solutions of 
the field equations in terms of one arbitrary function of £ (|32 ]) - ([34l) . 

These solutions are Riemann-flat in 5D but curved in AD (|44p . This latter property, together with the property 
(i) mentioned above, allow us to interpret these cosmologies as Riemann-flat embeddings for spatially-flat FRW 
cosmologies in AD. In this interpretation our universe arises as a spontaneous separation from an empty 51? Minkowski 
space, as envisioned by Zeldovich [7]. 

7 This is equivalent to requiring continuity of the first and second fundamental forms across a 3D surface r = tq = constant 1451 . 
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We have analyzed in some detail the AD interpretation of ([52")) (for 1 = 2), which admits the simply homothethic 
Killing vector £ A = (t, r, 0, 0, y). We have seen that the AD projection of this vector does not constitute a 
AD nomothetic Killing vector for the metric induced on Y, y (|43p . For the physics in AD this opens a wide range 
of non-homothetic, but still self-similar, possibilities for the evolution of our universe (|A-ip - i|A-3|) . among them it 
allows the introduction of a cosmological constant and a de Sitter phase of exponential expansion. 

The question is how to establish the arbitrary function of £ in (|32"f - (|34p . There are many plausible ways for 
selecting it. As an example one can choose some "geometrical" criteria, e.g., choose a particular coordinate frame or 
to assume some specific symmetry, as we do in sections 4.1, 4.2 and the Appendix, respectively. A distinct approach, 
which is probably more satisfactory from a "physical" point of view, is to select the arbitrary function by imposing 
conditions on the effective picture in AD; this is what we do in section 4.3. Certainly, the existence of an arbitrary 
function allows us to accommodate a number of physical models. 

The model discussed in 4.3 roughly presents the essential features of the evolution of our universe. Specifically, it 
shows inflation from an original birth, not described by the de Sitter solution (|7D|) , followed first by a FRW phase and 
then by an accelerated expansion. The effective cosmological constants during inflation A( in ^) and after inflation A 
are not equal, otherwise the universe cannot enter a phase of accelerated expansion (|76p . 

In summary, in this work we have studied a family of self-similar cosmological metrics in 5D. As far as the author 
knows, this is the first work where self-similarity is used for finding exact solutions to the field equations in 5D. 
The solutions may be used as 5D embeddings for the spatially-flat FRW cosmological models of ordinary general 
relativity in AD. 

We have not investigated fully the possible physical interpretations of (f3"2"| . we have just considered a particular 
case. Neither have we investigated the solution (f34|) . which seems to be totally different from (|32|) . An immediate 
extension of this work is the study of self-similar cosmologies with extra dimensions without the assumption of spatial 
isotropy and/or homogeneity. Also, it is interesting to find the transformation of coordinates from (|3"2")) - (f3~4"|) to 
given by (|36|) . This is a non-trivial task, but the existence of such a transformation is guaranteed by the fact that 
the Riemann tensor in 5D vanishes. With the appropriate transformation of coordinates at hand one can extend 
and generalize previous investigations [28] and study the birth of the universe in more detail. 



Appendix: Homothetic symmetry on E z 



Our aim here is (i) to show that the requirement of homothetic symmetry on E y singles out one specific metric in 
AD, namely, the spacetime part of the 5D metric ([5]), and that (ii) the homothetic vector is not parallel to (|43|) . 
In order to do this, let us express the induced metric on T, y , which is given by (|42"|) , in terms of S = e A / 2 , 



ds 2 = h^dx^dx" = ^fdt 2 - S 2 (£) [dr 2 

The Lie derivative of this metric along the 4D vector 

C = (At, Br, 0, 0), 

where A and B are some constants, is given by 



r 2 dfl 2 ] 



(A-l) 



(A-2) 



Thus, the requirement 

generates two independent equations, viz. 



£-c,hij = 2hi 



C^hap — 2h a p, 



A 



S 



B 



(A-3) 
(A-4) 



_ 2(1 -A) flnH S 6 (1-1?) 

57 -~az~' ^-^r- 



(A-5) 
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From the first equation we get S = C\£} 2 A ^ A + C2, while from the second one we obtain S = Cs^ 1 B ^I A , where 
Ci, C2 and C3 are constants of integration. Compatibility of these expressions demand (A — B) = 1, and C2 = 0. 
Consequently, S ~ ^ 2 - A )/ A and e" ~ £2(i-A)/A 

Since A is an arbitrary parameter, without loss of generality we can set 

2a- 1 , . . 

A= . (A-6) 

a 

With this selection the metric in AD becomes 

ds 2 = A e {l ~ a)/{2a - 1] dt 2 - Be /{2a - 1] [dr 2 + rW] , (A-7) 
which is identical to the spacetime part of © and admits the homothetic Killing vector 

C=f^*,^r,0,o). (A-8) 
\ a a J 

Clearly, this vector is not parallel to £g = (t, r, 0, 0). We emphasize that the 5D metric Q is homothetic along the 
5.D vector £ A = (t, r, 0, 0, y), but its spacetime part is homothetic along (|A-8[) and not along ^ = (t, r, 0, 0). 
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